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We reconsider the possibility that the observed baryon asymmetry was generated by the 
evaporation of primordial black holes that dominated the early universe. We present a 
simple derivation showing that the baryon asymmetry is insensitive to the initial black hole 
density and the cosmological model but is sensitive to the temperature-dependence of the CP 
and baryon- violating (or lepton- violating) interactions. We also consider the possibility that 
black holes stop evaporating and form Planck-mass remnants that act as dark matter. We 
show that primordial black holes cannot simultaneously account for both the observed baryon 
asymmetry and the (remnant) dark matter density unless the magnitude of CP violation is 
much greater than expected from most particle physics models. Finally, we apply these 
results to ekpyrotic/cyclic models, in which primordial black holes may form when branes 
collide. We find that obtaining the observed baryon asymmetry is compatible with the other 
known constraints on parameters. 
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1 Introduction 



Evidence for matter-anti-matter asymmetry is inferred from studies of the particle con- 
tent in cosmic rays, the absence of an intense 7-ray background from matter-anti-matter 
annihilations on cluster scales pQ, big-bang nucleosynthesis (BBN) [2 J and high precision 
measurements of the cosmic microwave background (CMB) [3J. However, imposing the ob- 
served baryon-to-entropy value, B — ^ = nb ~ nb = 9. 21^4 x 10~ n , as an initial condition 
seems extremely unnatural. The currently favored alternative is to imagine that the baryon 
asymmetry was generated from symmetric initial conditions through a dynamical mecha- 
nism. 

The idea of producing an excess of baryons over anti-baryons in black hole evaporation 
goes back to Hawking [4] and Zeldovich [5J. More recently, this notion has been explored 
extensively in the context of grand unified theories (GUTs) (e.g. [Hid IE])- As a black hole 
evaporates, its horizon temperature eventually becomes higher than the GUT-scale. Hawking 
evaporation of a gas of primordial black holes provides out-of-equilibrium conditions, whereas 
grand unified theories generically violate CP and baryon number conservation. In particular, 
evaporating black holes may emit massive particles whose decay generates a net baryon 
asymmetry. Hence, black hole evaporation and local GUT-scale physics together satisfy the 
Sakharov conditions [9] and may provide a natural mechanism for generating the observed 
baryon asymmetry of the universe. Alternatively, it is possible to envisage generating a net 
lepton number L at lower energies and subsequently producing the baryon asymmetry via 
B-L conservation [TIH[TT1[T2]. 

In this paper we reconsider direct baryogenesis or baryogenesis via leptogenesis from the 
evaporation of primordial black holes. Under the assumption that black holes dominate the 
energy density of the universe by the time they evaporate, we present a simple derivation of 
the resultant asymmetry that is insensitive to the cosmological model. Hence, we can apply 
our results both to inflation [T3] and ekpyrotic/cyclic models [HI [T5]. We find that the 
baryon asymmetry is independent of the initial black hole number density, but does depend 
sensitively on the nature of the CP and B or L-violating interactions. We also consider the 
possibility that black hole evaporation stops due to quantum gravity effects when the mass 
reaches the Planck scale and the black holes form massive relics that act as dark matter. 
The cosmological consequences of Planck- mass black hole relics have been studied previously 
by MacGibbon [16J, Carr et al. [T7] and Barrow et al. [18]. Here we discuss the combined 
constraints if black holes form relics and also account for baryogenesis. 

What processes in the early universe could have generated primordial black holes? Pri- 
mordial black hole formation from big bang inhomogeneities was first discussed by Carr and 
Hawking [T3]. However, if inflation removes all pre-existing inhomogeneities, any cosmologi- 
cally interesting black hole density has to be created after inflation. A number of mechanisms 
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have been proposed for generating post-inflationary primordial black holes. Carr [20J dis- 
cussed black hole formation if the spectrum of density fluctuations generated during inflation 
is very blue (spectral index n s > 1). If the primordial power spectrum is nearly-scale invari- 
ant and doesn't have significant running, then this mechanism is now tightly constrained 
by CMB measurements [3]. Primordial black hole production has also been proposed in 
the context of hybrid inflation. Garcia-Bellido et al. [21J showed that density perturbations 
in hybrid inflation models can reach a large enough magnitude to produce black holes for 
wave numbers of order the horizon scale during the transition between the two inflationary 
stages. The mass and number density of black holes produced in these models is highly 
model-dependent, including large regimes which are cosmo logically disastrous or insignifi- 
cant. Hsu [22] studied the formation of black holes in models of extended or hyperextended 
inflation (see [23]). These models contain two sources of fluctuations - the usual nearly 
scale-invariant scalar field fluctuations and also inhomogeneities created by the collision of 
bubbles |24j produced during the last few e-foldings. Barrow et al. [7] provided an analytical 
treatment of GUT baryogenesis from these black holes. 

In this paper we present a model-independent derivation of the baryon asymmetry gener- 
ated by primordial black holes that can be applied in both inflationary and ekpyrotic/cyclic 
models. In the cyclic scenario black holes may form when branes collide. Since the reheat 
temperature in ekpyrotic/cyclic models is generally modest, it is important to have a re- 
liable baryogenesis mechanism operating at temperatures well below the GUT-scale. We 
show that black hole baryogenesis provides a mechanism that is effectively independent of 
the reheat temperature of the universe. The potential challenge is that the ekpyrotic/cyclic 
model parameters must also satisfy a number of other conditions in order to resolve the 
horizon, flatness and monopole problems and to obtain an acceptable spectrum of density 
perturbations [25]. However, here we show that black hole baryogenesis is compatible with 
all current constraints for a wide range of parameters. 

In both inflationary and ekpyrotic/cyclic models, precise calculations of the initial black 
hole density are either difficult or highly dependent on model parameters. One of the im- 
portant features of the scenario we are considering is that the baryon and possible remnant 
densities are insensitive to the initial black hole density provided the black holes dominate the 
universe by the time they decay. Reliable predictions for the baryon asymmetry can there- 
fore be made even in the absence of an exact calculation of primordial black hole production. 

The outline of the paper is as follows. In section [21 we examine black hole baryogenesis. 
After considering Hawking emission from a single black hole, we compute the baryon number 
generated by an ensemble of black holes. We explore the parameter space of CP and B- 
violating interactions that reproduce observations. Section [3J applies our results to black hole 
baryogenesis in the ekpyrotic/cyclic model. In section 0] we briefly discuss the cosmological 
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implications if Planck mass remnants form as the final stage of black hole evaporation. We 
summarize our conclusions in section 

We use natural units, h = c = k B = 1, throughout, and define SnG = l/m P1 = 1. 
For numerical estimates we use the reduced Planck mass mpi ~ 10 18 GeV and Planck time 
tp\ = 1/mpi k5x 10~ 44 sec. 



2 Black Hole Baryogenesis 

We begin by giving a brief review of baryogenesis during Hawking evaporation [26] of a single 
black hole (BH). Then we consider an ensemble of black holes and compute the resulting 
baryon asymmetry. 

The horizon temperature of an uncharged, non-rotating Schwarzschild black hole is 

1 



BH 



M BH 

Its rate of evaporation is given by l27j.Fl 



dM BU _ vr 2 4 7T g 

~dt ^* BH bh " -480 m ? 



9*AbhT bh - -j^T7T~ ' ( 2 ) 

BH 



where A B ^ — 47ri? BH is the surface area of the black hole event horizon, M B n(t) is the 
instantaneous black hole mass and g* ~ 100 is the effective number of relativistic degrees of 
freedom of particle species radiated at temperature Tbh- Integrating equation (j2J) gives 

M BH (t) = M n--J , (3) 

where the black hole lifetime is 

160 1 _ _o . iS 

T= M 3 , (4) 

7T g* 

for a black hole with initial mass M . The Gibbons-Hawking temperature (JTJ) implies the 
differential mass decrease 

dM BH = --^dT BU = -dE , (5) 

BH 

where dE is the energy emitted by the black hole when its mass decreases from Mbh to 
M B h + dM BB . Since the radiated particles have mean energy 3T B h, the differential number 
of particles emitted is 

dE 1 1 

dN=-— = --^dT BR . (6) 
<jJbh <j j bh 



^^Note that Barrow et al. [7] use an approximate equation for the mass loss of a black hole which under- 
estimates the true rate by a factor of 10 5 . This leads to an overestimation of the final baryon asymmetry by 
a factor of 100. 
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Integrating equation (jH]) beginning from the initial black hole temperature T = Tbh (Mo) = 
1/Mq we find the total number of particles emitted by a single black hole of initial mass M 

poo jif2 

N= / dN(T) = -±. (7) 

Baryogenesis can occur through the emission of a particle species X with mass Mx, whose 
subsequent decay violates baryon number conservation. If the initial black hole mass is small 
enough, T = 1/M > Mx, the fraction of X particles emitted is given by the equipartition 
factor 

fx~—, (8) 

9* 

where gx denotes the number of degrees of freedom of particle X. The total number of 
.B-violating particles is 

N x = f x N if M <-^. (9) 

Mx 

Significant X-emission requires Tbh > Mx- Very massive black holes are initially not hot 
enough to emit X particles (T < Mx), and baryon production during the initial stages of 
evaporation is exponentially suppressed. In this case, the lower limit (T ) in the integral in 
equation (J7|) should be replaced by Tbh = Mx, and equation ([9]) becomes 

N X = f x r dN(T) = f x ( 1 -^-) 2 ^- if M >-L. (10) 



Hence, we have 
where 



Mx \ M o J 6 M x 



N x = Kf x N, (11) 



, 1 HM <1/M X 
~ i {M M X )- 2 if M > 1/Mx ' 



.B-violating X-decays may still lead to a baryon symmetric universe if the corresponding 
anti-particle processes occur at the same rate, i. e. baryon asymmetry requires CP violation. 
If the CP violating parameter is 

_ ^ T(X -> ft - T(X -> j) 
7 = 2_.Bi , (13) 

where Bi is the baryon number of the final state fa, then the excess baryon number generated 
by the evaporation of a single black hole is 

jN x = jKfxN . (14) 
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Next, we consider the baryon asymmetry generated by an ensemble of small black holes. 
Imagine a dilute gas of primordial black holes, formed through some mechanism that we 
leave unspecified for the purpose of this general discussion. Depending on the details of 
the cosmological model, the universe may also contain a background of radiation and other 
contributions to the total energy density. As we shall see, the resulting baryon asymmetry is 
independent of these details provided that the black holes survive long enough to dominate 
the universe by the time they evaporate. 

In some of the cosmological models of interest, including the extended inflation and 
ekpyrotic/cyclic models, black holes are produced with a narrow distribution of masses set 
by the horizon size at the time of formation. Consequently, we assume for simplicity an 
initial distribution in which all black holes have the same mass M . (It is straightforward to 
generalize our results to a broader distribution.) Equation (T5]) shows that the most energy 
is transfered as the mass of the black hole becomes small during the final stages of the 
evaporation process. This and the fact that ^formation *C r (where r is the mean black hole 
lifetime) allows us to treat the evaporation as if all particles were produced at a single 
instant, t ev ap ~ T - Approximating equation Q by a step-function, M B h(£) ~ M 6(r — t), 
and assuming that black holes dominate the universe at the time they decay, we find 

Pevap = PBH(Cap) ~ M , (15) 

where = nBu(t cvap ) is the black hole number density at t evap . If the black holes dominate 
at evaporation then the total entropy and baryon density of the universe is determined by the 
black hole evaporation products. All black hole energy density is transformed into radiation 
at tevap • Hence, we have 

PBH(Cap) = Ay(^ap) ■ ( 16 ) 

The radiated particles equilibrate with the surroundings and the temperature of the universe 
is defined via p^ vap ~ T e ^ ap . Using equations (IT51) and (ITtjj) we have 

<n P M = ^T e 4 vap • (17) 

The Hubble parameter at evaporation (black hole domination) is given approximately by 

2 1 

#evap ~ 7TT > (18) 

" ''evap 

where t evap is the FRW time. The last approximation is valid provided that t ev ap ^ 
^BH-domination- Using the black hole lifetime (j4]) and H^ vap = p C va P /3 in combination with 
equation (flBT) gives 

#2 = 1 (JL) 2 £ = - n e Z P M . (19) 
evap 9 V160/ M n 6 3 K J 
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Solving equation ffT^l) for the black hole number density at evaporation we get 

Note that the black hole number density at evaporation is independent of both the number 
density at formation and the cosmological model. This is a direct consequence of the as- 
sumption that the black holes come to dominate the universe at some point in the evolution. 
It means that the cosmological conditions at evaporation that determine both the entropy 
and the baryon number are characterized by a single parameter, M . 

From equation ( jTTl) we then obtain the evaporation temperature 

T = (l.) 1 ".(Jt-Y /i .-L- (21) 
™* \32j UOoJ Mf-' 

In order to reproduce the success of standard big-bang cosmology, we require the universe to 
be radiation dominated at nucleosynthesis. The black holes therefore must have evaporated 
by that time, which leads to the following BBN constraint 

T eva ,p > T B bn ~ 10~ 22 . (22) 
As a limit on the initial black hole mass, this gives 

M <4xl0 14 (^) 1/6 . (23) 

Using equation (fT4l) for the net baryon number created by a single black hole, the number 
density of baryons at the time of evaporation is 

n B = n B (t evap ) = jKfxNn^ ■ (24) 

We define the effective entropy radiated by a single black hole as 

_ s 

°~ = evap ) (,25) 
n BH 

where the entropy density generated by a gas of black holes is 

x 2tt 2 o 

s = s(t cvap ) = —g*T cvap . (26) 
Combining ff24|) and (T26]) . the baryon-to-entropy ratio is 

B = ^L = mMi. (27) 

s o 
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Equations (J2"5j) . (J2Uj) and (I2T1) give the entropy produced per black hole 



3T cvap 3\5j UoO/ 1 ' 

This shows that the effective entropy radiated by a single black hole is a factor of \[Mq 
larger than the initial Bekenstein-Hawking entropy of the black hole, S = ^/^^4bh ~ Mq- 
Substituting N = Mq/Q and ct(Mq) into (1271) . we finally get our master expression for the 
baryon number of the universe 



B 



uoo/ 



k . (29) 



M 1/2 



The parameters g*, 7, /x, and Mx are set by microphysics, whereas M is determined by the 
cosmological model. Using f x ~ l/g* the first factor in (T2"9"j) is C(<7*) = 1CT 3 (^) 3//4 . Note 
that B is independent of the initial black hole number density. This illustrates that the final 
baryon asymmetry is insensitive to the initial conditions as long as evaporation occurs after 
the black holes dominate. The final baryon asymmetry (T29j) displays a simple dependence 
on the free parameters Mo, 7 and Mx- B scales linearly with 7, the net baryon number 
created by each black hole. For small initial black hole masses, M < 1/Mx, k — 1 and 

1 /2 

B is simply proportional to the ratio 7/M . Increasing M therefore suppresses baryon 
production. Very large initial masses, M > 1/Mx, lead to a further suppression from the 
factor k = 1/(M M X ) 2 < 1. 



Figure[T]is a plot of [5=10~ 10 ]-contours in the 7-Mx plane for three different fixed values 
of the initial black hole mass Mq. For Mx < 1/Md (immediate emission of X particles) 
7 = ^Mq 1 ^ 2 = const, independent of Mx, giving a horizontal section. For more massive X 
particles, M x > 1/M , baryogenesis is suppressed and 7 has to be increased to compensate 
that suppression. This explains the rise of the contours for large Mx- The transition between 
these two regimes occurs at Mx = 1/Mo- Each contour corresponds to combinations of the 
three parameters 7, Mx and Mq that reproduce the observed baryon asymmetry, B pa 10~ 10 . 
The region below the contours corresponds to parameter combinations that are unacceptable 
because they underproduce baryons, B < 10 -10 . This is the forbidden region. The region 
above the contours overproduces baryons, B > 1CT 10 . We call this the allowed region. 
Overproduction is allowed since we can imagine non-adiabatic processes after the black holes 
evaporate that create extra entropy and dilute the abundance to the correct value. We see 
that there is a wide range of allowed parameters. We also notice that increasing Mq decreases 
the allowed region. Small masses, therefore, seem to be preferable, although consistency with 
the assumption that the black holes live long enough to dominate the universe at evaporation 
has to be tested separately. This can only be done once the cosmological setup is specified. 
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Figure 1: [5=10~ 10 ]-contours in the j-Mx plane for three different values of the initial black 
hole mass M (cf. equation (1251) ). 7 determines the amount of CP violation, while Mx is the 
mass of the B- violating particle X. The contours correspond to parameter combinations that 
produce B pa 10~ 10 . Regions below the contours lead to a baryon number that is too small 
to be consistent with observations (B < 1CT 10 ). Regions above the contours overproduce 
baryons (B > 10~ 10 ). Since entropy production at a later stage of evolution may dilute B, 
this region is viewed as a cosmologically acceptable region. The dotted line corresponds to 
the leptogenesis model of [28] which predicts 7 = 10 2 Mx- The dashed line defines the region 
of parameter space which, in principle, allows black holes to simultaneously be the source of 
baryons and dark matter (see section H]). Note that this requires very large CP violation for 
small Mx- 

Thus far, we treated 7 and Mx as independent free parameters. We now show how this 
plot can be used to explore a given particle physics model of baryogenesis that relates 7 and 
Mx- For example, we consider the leptogenesis model in [28] that relates the CP violation 
parameter 7 to M X - 

7 = 10 2 M X - (30) 

Using the above analysis we can easily evaluate the feasibility of this black hole leptogenesis 
scenario. The dotted line in Figure [H corresponds to equation (1301) . For small M there 
is a finite range of j(Mx) in the region that allows sufficient baryogenesis. In Figure [H 
this is illustrated for M Q = 10 4 . The dotted line crosses the [B = 10~ 10 ]-contour at two 
points. These points correspond to combinations of j(Mx) and M Q values that reproduce 
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observations precisely. In between these two points there is a finite line segment where 
baryons are overproduced, which is also allowed assuming some modest dilution later. For 
this particular black hole leptogenesis model to be feasible, we conclude that 

M < 10 6 ; (31) 

otherwise, the dotted line will lie completely in the forbidden region and baryons are under- 
produced for all choices of parameters. 

3 Black Hole Baryogenesis in the 
Ekpyrotic/Cyclic Models 

Our derivation of black hole baryogenesis is insensitive to the cosmological model provided 
that the black holes dominate the energy density of the universe at the time of evaporation. 
An interesting application is to the recently proposed ekpyrotic/cyclic models [IH[T5]. Ac- 
cording to the ekpyrotic model the hot big bang is caused by the collision of two 3-branes 
bounding an extra spatial dimension [H]. In the cyclic scenario this may repeat itself peri- 
odically, leading to a periodic sequence of big bangs and big crunches. The model is designed 
to resolve the cosmological horizon, flatness and monopole problems and generate a nearly 
scale-invariant spectrum of density perturbations without invoking a period of high-energy 
inflation [13]. The ekpyrotic/cyclic model may also provide a natural setting for black hole 
baryogenesis. The ekpyrotic mechanism is designed to smooth the universe on superhorizon 
scales before the bounce occurs. At present, a bounce is not proven to occur, but, assuming 
it does, it is reasonable to suppose that instabilities develop within the horizon and produce 
black holes on subhorizon scales. In this section we consider whether black hole baryogene- 
sis is compatible with the constraints on ekpyrotic/cyclic models that are already known [25J. 

According to the cyclic model, the universe is dominated by the brane kinetic energy 
immediately after the bounce. In the effective 4D theory, the separation between the branes 
is described by a scalar field <fi and the brane kinetic energy is characterized by the scalar 
field kinetic energy oc a(t)~ 6 . A generic potential for is shown in Figure [2j Region 
I corresponds to the present epoch, which is characterized by dark energy domination and 
accelerated expansion. Scale-invariant density perturbations are created during a phase of 
slow contraction (region II). Chaotic behavior in the big crunch is suppressed by the expo- 
nential form of the potential (w>l, see [29]) and fluctuations are linear until t^ d . However, 
short- wavelength fluctuations become highly non-linear after t^ d (region III). The distance 
between the branes becomes small and quantum gravity corrections become important. Clas- 
sical general relativity is insufficient to describe the further evolution. It is known, however, 
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Figure 2: Black hole baryogenesis in the ekpyrotic/cyclic universe. The scalar field pa- 
rameterizes the distance between the branes. V(<p) is a generic potential for 0. Region I 
corresponds to the observed dark energy domination today. A scale-invariant spectrum of 
perturbations is generated in region II. Perturbations are stabilized until t^J, after which 
local instabilities develop and form black holes. The size of these black holes is limited by 
the horizon size at t[ n \. The cosmological implications of these black hole for the post-big 
bang evolution of the universe are studied in this paper. 



that the effective gravitational constant becomes large when the branes collide and particle 
collisions reach large energies. The horizon at time t sets the maximal distance a signal sent 
at time t can travel prior to the bounce. We therefore assume that the instabilities that 
develop in the big crunch are limited by causality to small scales, the maximum scale being 
the horizon at t^} (the time when perturbative control is lost). This sets the maximum 
black hole mass, but we don't call it a black hole before the bounce, since the cosmologi- 
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cal horizon shrinks to become smaller than the Schwarzschild horizon. The horizon grows 
after the bounce to encompass the instability at about t^d ( see Figure [2]). We define black 
hole formation as the time tf after the bounce when _f/" _1 (£f) = H^ d = H -1 ^^), *.e., the 
horizon becomes comparable to the Schwarzschild radius. The near symmetry of the time 
evolution of the scalar field (H(tl~l) « H(t^l)) then allows us to use the horizon at t[~l as 
a scale for the black hole formed at tf ~ t^\. 



log(p) 




If t r t bh tevap log(a) 



Figure 3: In the cyclic model, the universe changes from scalar field dominated (0D) when 
the black holes form (at t — tf ~ ^nd)' ^° radiation dominated (RD) at t r , to black hole 
dominated (BHD) at t hh , and back to radiation dominated at t C va P - 

Here we focus on the evolution after the bounce. Assuming the collision produces signif- 
icant amounts of radiation and black holes, the various densities evolve according to Figure 
[3j The case we study in detail below (and that is illustrated in Fig. [3]) assumes for simplicity 
that the initial radiation density produced at the collision overtakes the scalar field kinetic 
energy before the black holes. This is equivalent to assuming that the energy density in black 
holes at radiation-scalar field equality (t r ) is less than the energy density in radiation. We 
later reformulate this as a constraint on the number of black holes per horizon at the time of 
formation. Since radiation dilutes more quickly than the pressureless gas of black holes, the 
black holes will eventually come to dominate the universe (t > tbh)- As we have shown, the 
baryon asymmetry resulting from their evaporation is then independent of the initial black 
hole number density, but depends sensitively on the nature of the CP and B or L-violating 
interactions. The initial radiation energy density redshifts away in the subsequent black hole 
dominated phase and plays no significant role afterwards. Black holes dominate for a finite 
time until they decay at t ev ap- Evaporation occurs before nucleosynthesis (t ev ap < ^bbn), so 
standard cosmology is recovered after that time. 

The previous discussion suggests that the bounce produces both a small amount of radia- 
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tion from the non-adiabatic collision of the branes and a gas of primordial black holes. Once 
the radiation or black hole dominated epochs begin, the scalar field kinetic energy is rapidly 
damped and the field comes to halt. In order to cycle, radiation or black hole domination 
should not begin until has had time to go from — oo, past V = — \ V en d\, across the potential 
well and back up the potential to the plateau where V = O(V ) (see Fig. [2]). Otherwise the 
scalar field would get trapped in the minimum at V en d and the universe would be anti-de 
Sitter space. Khoury et al. |25J showed that this leads to the following cycling constraint on 
the radiation energy density (or temperature) at scalar kinetic energy-radiation equality (t r ) 



Pr ~ T r 4 < Knd -M , (32) 

V ^cnd / 

where e is the fast-roll parameter defined in ekpyrotic cosmology [15J. Under the assumption 
that radiation from the collision dominates before the black holes, we have 

Per < P7 ~ Pr ■ (33) 

This can be re-expressed as a condition at the time of black hole formation. The black hole 
number density at radiation-scalar field equality is related to the density at formation via 

-bh = -bh(^) 3 = 4h^, (34) 

where we used H 2 oc a~ 6 during domination of scalar field kinetic energy (tf ~ — » t r ). 
Assuming that the initial black hole mass is equal to the horizon mass at formation, M ~ 
PfHj 3 ~ Hj 1 , and using H 2 ~ T r 4 gives 



V 

n 



bh - n f BH M T? . (35) 
Defining the number of black holes per horizon at formation as 

N B r = n f BU Ml , (36) 

we find 

^BH = A^BH " || • (37) 

Imposing the cycling constraint in the form pg H ~ ng H Mo < pL ~ T r 4 , we obtain 

iV B H < M T 2 , (38) 

where T r obeys the cycling condition ( 1321) . Black hole baryogenesis works for larger A^bh> but 
then black holes dominate before the collisional radiation and a separate analysis is needed. 
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We also impose the constraint that the black holes survive long enough to eventually 
dominate the universe. This is required for the results of section [2] to be applicable. We 
denote quantities evaluated at radiation-black hole equality by (bh). The condition that the 
black holes survive long enough to dominate the universe is r > t^, or in terms of the scale 
factor of the universe 

Oevap > «bh • (39) 

At radiation-black hole equality we have, psH^bh) = P7^bh) , or 

"BH^o ( = T r 4 ( ±X . (40) 

\ «bh / V «bh / 



From (j40p and (157)) we obtain 

Or_ = ^BH _ (41) 

a hh M T? 

This confirms that the cycling condition (13"51) is equivalent to abh > <V Using a r /a hh = T hh /T r 
we get the temperature of the universe at radiation-black hole equality 

T - = lt' < 42 > 

During black hole domination (tbh — » t cvap ) we have 

„ \ 3 tt2 rp4 at4 

6Vap ^ " bh " bh -Ml^-. (43) 



a hh J l/r» U M 4 T* 



evap 

Notice that -ff ev a P is related to the lifetime of the black holes. The condition that the black 
holes survive long enough to dominate the universe, a ev a P > ^bh, translates into H evSuP < H^, 
or 

iV BH > • (44) 

This gives a lower limit on the initial black hole mass 

M ">W, (45) 

where T r is bounded by the cycling condition. T r is determined by the initial ratio of 
scalar kinetic energy and radiation density, with larger T r corresponding to more radiation 
produced at the bounce. Keeping everything else fixed, while increasing T r , increases the 
time to radiation-black hole equality. As fj45|) shows, the black holes then have to be more 
massive, in order to survive long enough to dominate. A^bh is related to the number of 
black holes at formation. For a large number of primordial black holes, radiation-black hole 
equality occurs earlier. Hence, increasing A^h while keeping T r (or the initial amount of 
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radiation) fixed, decreases the time to radiation-black hole equality. The constraint on M 
therefore weakens. 

Combining the cycling (|38|) and survival (|44p conditions gives 

T, 



T — < N BU < M T? , (46) 



where 



Hence, cycling and survival can simultaneously be achieved only if the following consistency 
constraint is satisfied 

M > T r - 2/3 . (48) 

As discussed earlier, the maximal initial black hole mass is set by the horizon mass at the 
end of the ekpyrotic phase: 

M < M Hor (4<l) ~ Piinl) l^endP 3 ~ l^end^ • (49) 

In [25] it was shown that the horizon at t cn d is 

l^endP 1 ~ (2eV cnd )- 1/2 ■ (50) 
This gives the following relation between M , e, and V en <i 

M <(2eKnd)- 1/2 . (51) 
Hence, we have related the initial mass of black holes to the depth of the ekpyrotic potential. 



Figure H] (adapted from [25]) illustrates the region of parameter space for which the cyclic 
model satisfies all known cosmological constraints. The outer triangular region in the V en< i-T r 
plane encloses the range of parameters that satisfy the BBN, cycling and gravitational wave 
background constraints studied in [25]. For the purpose of illustration, assuming that the 
black hole masses saturate the bound ( 15T1) . i.e. M = (2eV en d)~ 1 ^ 2 , we consider the constraint 
set by black hole baryogenesis. We impose an upper bound on CP violation, 7 < 10 _1 . From 
the analysis presented in section [2] (7 m i n ~ 10~ 7 M^ 2 ) this implies M < 10 12 . Using the 
relation between M and V en( ± we get the baryogenesis limit 

> {2eY 1/i 1(T 6 w 1(T 6 = 10 12 GeV , (52) 

for e ~ 10~ 2 . In addition, we can re-express the consistency condition (1481) 

V^<(2e)-^Tl'\ (53) 
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Figure 4: Constraints on the cyclic model with (grey triangle) and without (outer triangle) 
black hole baryogenesis: The outer triangular region (reproduced from [23]) shows the range 
of V^ 4 and the reheat temperature T r , over which the cyclic model satisfies all known 
cosmological constraints, with fixed e = 10~ 2 . Relating the initial black hole mass Mq to the 
horizon mass at the end of the ekpyrotic phase leads to an estimate of Mq in terms of the 

— 1/2 

depth of the ekpyrotic potential Knd : M ~ V cnd . The black hole baryogenesis analysis of 
section [2] is only applicable if M > T^ 2 ^ 3 , or < Tr . This is labeled the consistency 
constraint. Imposing an upper bound on the CP violation parameter 7 < 10 _1 and using 
7min ~ 1CT 7 Mq^ 2 from section [2] gives the baryogenesis limit: M < 10 12 or V^ 4 > 10 12 
GeV. The grey region corresponds to parameter combinations for which primordial black 
hole baryogenesis is viable and consistent with known constraints on the cyclic model. 



Equations (|52|) and (|53|) together with the constraints from [25] define the grey region in 
Figure HI This corresponds to the region of parameter space for which the black hole baryo- 
genesis mechanism of section [2] is viable. We point out that regions outside the grey area are 
not strictly disallowed, but correspond to regions of parameter space that either require a 
separate analysis or where insufficient information is available to assess the viability of black 
hole baryogenesis. The consistency constraint is not a real physical constraint, but merely 
the limit of applicability of our analysis. In that sense the shaded area may be viewed as 
a minimal allowed region. Finally, the allowed region should be interpreted in combination 
with the generic results of section [2] (Figure [Q. These imply that large M (small Knd) 
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may require unconventional particle physics to give the large value of 7 required to obtain 
B = 10 -10 . Successful black hole baryogenesis is therefore more physically plausible in this 
sense towards the left (increasing V en< x /decreasing M ) of the allowed region. 



4 Black Hole Remnants as Dark Matter? 

Hawking's semi-classical treatment of black hole evaporation breaks down as the mass of the 
black hole approaches the Planck mass and quantum gravitational effects become large. A 
complete theory of quantum gravity seems necessary to decide on the final state of black 
hole evaporation. So far we have assumed that the black holes evaporate completely and 
the final black hole temperature tends to infinity. We now extend the discussion to include 
the possibility that the black holes cease evaporating when the mass is of the order the 
Planck scale due to quantum gravity effects and form stable Planck mass relics. Whether 
this concept is plausible or not is beyond the scope of this paper (but see [35]). Here we 
simply consider whether such remnants, if they form, could account for the dark matter. 



If fij is the ratio of the density of component i to the critical density, then dark matter- 
radiation equality occurs when = ° M o ° ^ = 1, and the temperature at equality is 

rp _ ^DM,0 rp (r A \ 

-tDM=7 — -pj J-0 , \p^) 

"7,0 

where To is the CMB temperature today. From black hole evaporation, the densities of 
radiation and relics are 

P7el P = ^BH PM rel, (55) 
p evap = n evap Mo _ (5g) 

Note that in (1551) we assume for simplicity that the initial mass spectrum has a very narrow 
range about Mq and can therefore be approximated by a ^-function. 

evap 

Equations (1551) and (I5BT) imply relic-radiation equality at %kp arcl=7 = 1, with corresponding 
temperature 

T re i=7 , ~r T cva p . (57) 

Mq 

If the black hole relics are the dark matter, then 

T rel=7 = T DM=7 , (58) 

or, from equations f!5^|) and fl57j) . 

^DM,0 _ ^rel T evap 



Oy,0 M T 
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Substituting T evap from equation (J2"lj) gives 

fi 7 , To \32j 1 100/ m 5/2 ' 1 ' 

Assuming M rel ~ 1, o* ~ 100, and the WMAP results G2, « 10 4 and T = 2.7 K ~ 

10~ 31 , we obtain 

(M ) DM ~ 10 11 . (61) 

Black hole relics can be the dark matter if the characteristic initial black hole mass has 
this virtually unique value. If M < 10 11 then fl ve \ > Odm, and relics are overproduced 
compared to what observations allow. Hence, if stable black hole relics are the final stage 
of Hawking evaporation then their initial mass must be greater than (Mo)dm = 10 11 , an 
important constraint on black hole baryogenesis. 




Figure 5: Constraints on the ekpyrotic/cyclic model: Same as Figure HJ but with relic dark 
matter limit included. The grey region is the new allowed region. 

It is interesting to consider whether the black holes could simultaneously be the source 
of dark matter and baryon asymmetry. For a generic black hole baryogenesis scenario this 
is easily analyzed by reconsidering Figure [H The dashed line in Fig. [Tjis the limit on 7-Mx 
for M = 10 11 . The region above the dashed line in principle allows primordial black holes 
to be the source of both the baryon asymmetry and the dark matter of the universe. Note, 
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however, that this requires a CP violation parameter of 7 > 10~ 2 for Mx < 1/Mq = 10~ n . 
Typical particle physics models of CP violation do not predict such large values for 7, which 
suggests that black hole baryogenesis and relic dark matter can't be realized simultaneously. 
Incorporating the analysis of black hole relics into our constraints on black hole baryogenesis 
in the cyclic model leads to Figure [51 Only the parameter space corresponding to large black 
hole masses and large 7 is now allowed. We emphasize that it seems very hard to conceive 
a realistic particle physics model of CP violation in this regime. 

5 Discussion 

Hawking evaporation of primordial black holes provides an interesting mechanism for gen- 
erating the baryon asymmetry of the universe. In this paper we presented a discussion of 
primordial black hole baryogenesis that is insensitive to the background cosmology and in- 
vestigated the parameter space for which the theory reproduces observations. A gas of small 
primordial black holes may be produced on subhorizon scales in the ekpyrotic/cyclic models. 
Applying our model- independent results for black hole baryogenesis to the ekpyrotic/cyclic 
model we showed that, for a large range of model parameters consistent with all other known 
constraints, these black holes could easily account for the observed baryon asymmetry. 

We comment briefly on the potential problem that any baryon excess generated prior to 
the electroweak era may be erased due sphaleron transitions [10]. The simplest possibility 
of avoiding sphaleron washout is to require that the black holes survive until after the 
electroweak phase transition, i.e. the reheat temperature is less than the electroweak scale, 
^evap < ?ew ~ 1CT 16 . This would imply a lower limit on the initial black hole mass, 

1 /6 

M > Ax 10 10 (^) 1 . This also forces us to a regime where 7 is very large. Because 
this corresponds to a relatively small range of parameters, we have considered instead the 
alternative of preserving the baryon asymmetry produced during black hole evaporation via 
B — L conservation [12]. This opens up a much wider and more physically plausible range 
of masses and couplings. 

Assessing the final state of black hole evaporation requires a better understanding of 
quantum gravity. If black holes evaporate completely, then their final temperatures become 
arbitrarily large. One might therefore worry about monopole production during the final 
stages of black hole evaporation. This was studied in [2H E21 E2], which conclude that 
this process is exponentially suppressed by semi-classical effects. We also investigated the 
possibility that black holes form stable Planck mass relics that could act as dark matter. 
We found that this would be very constraining for black hole baryogenesis. To avoid relic 
overproduction requires large initial black hole masses M > 6 x 10 10 . This, however, 
suppresses the final baryon asymmetry. Successful black hole baryogenesis therefore requires 
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relatively large CP violation 7 for very small Mx ~ 1/Mq. In principle, it remains conceivable 
to achieve simultaneously successful baryogenesis and a black hole relic abundance equal to 
the observed dark matter density, but only if CP violation parameter is much larger than 
most particle physics models predict. The more likely possibility according to our calculations 
is that relic production and black hole baryogenesis are incompatible. Either black holes are 
responsible for baryogenesis, or they form stable relics which could be the dark matter. 
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Note added: After completion of this paper, we received a preprint by Alexander and 
Meszaros [36] who considered similar issues but based their analysis on a different black hole 
baryogenesis mechanism, first proposed by Nagatani [37]. Nagatani considered generating 
the observed net baryon asymmetry through electroweak baryogenesis, in which a black hole 
with Hawking temperature above the electroweak scale creates a nearly static domain wall 
around itself through which a hot plasma flows; then, sphaleron processes within the domain 
wall generate the baryon asymmetry. Alexander and Meszaros present the optimistic view 
that this mechanism can produce the observed baryon asymmetry and simultaneously the 
observed dark matter density (assuming Planck mass black hole relics). However, Nagatani 
finds that the baryon asymmetry is suppressed by a factor of a^, where aw is the electroweak 
coupling, and the observed asymmetry can only be obtained if CP violation is maximal, 0(1), 
within the domain wall. This conclusion is similar to the case of the GUT baryogenesis 
mechanism we consider, and so does not change our qualitative interpretation that the 
observed baryon asymmetry and dark matter density are difficult to obtain simultaneously 
by the black hole evaporation mechanism. 
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